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Resumen 

En el presente trabajo se estudian dos metodologías distintas para la resolución de 
ecuaciones diferenciales parciales fraccionarias. La primera metodología consiste en 
usar la transformada de Laplace. Para esto, se propone una solución de tipo producto 
de funciones la cual nos permite reducir el problema a la resolución de dos ecuaciones 
diferenciales ordinarias fraccionarias y mediante el uso de la transformada de Laplace se 
obtiene una solución para las ecuaciones diferenciales obtenidas a partir del producto 
de funciones. Por otro lado, en la segunda metodología se asume que nuestras funciones 
admiten derivadas sumables y se realiza integración directa sobre la ecuacion diferencial 
parcial fraccionaria. Posteriormente se toma la transformada de Laplace del resultado 
obtenido y se llega a una ecuación de tipo Volterra de segundo orden. Finalmente, 
dada la convergencia de las integrales y funciones se simplifica la solución y se toma 
la transformada inversa del resultado previamente obtenido. Para finalizar el trabajo, se 
realiza una discusion sobre las limitaciones y libertades de ambos métodos. 


Palabras clave: Derivada Fraccionaria de Riemann-Liouville, Transformada de Laplace, 
Ecuaciones Diferenciales Fraccionarias, Ecuaciones Integrales, Derivada sumable, 
Ecuación integral de Volterra. 


Abstract 

The present work presents two methodologies to find solutions to fractional partial 
differential equations. The first methodology is an approach using the Laplace transform. 
In orderto do that, a solution as a product offunctions is proposed in order to separate the 
fractional partial differential equation in two fractional ordinary differential equations. Then 
the solutions of the fractional ordinary differential equations using the Laplace transform 


Editado por / are found. In the second methodology, we suppose that our functions admits a fractional 
a summable derivative in order to integrate them. Then, using the Laplace transform, a 
Blanco Volterra integral equation of second order is found. Using the convergence of the integral 
Recibido and functions the solution is simplified in order to take the inverse Laplace transform of it. 
Received: Finally, a discussion about the freedoms and restrictions oftwo methods is made. 
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INTRODUCCIÓN 


El modelado matemático consiste en traducir problemas relevantes de la industria y 
del día a día a problemas matemáticos mediante procedimientos teóricos o numéricos 
1]. En particular, el modelado de fenómenos físicos y químicos muchas veces nos lleva 
al planteamiento de ecuaciones diferenciales que pueden ser descritos mediante una 
función de una o más variables [2]. La idea del modelado es poder predecir y analizar el 
comportamiento del sistema; sin embargo, el modelado no describe en su plenitud el 
fenómeno que deseamos estudiar y una mejor descripción del fenómeno involucra un 
incremento en la complejidad del modelo [3]. Por ejemplo, consideremos un sistema 
masa resorte en el vacío, también conocido como oscilador armónico simple. A pesar 
de que el oscilador armónico nos ayuda a describir una gran variedad de sistemas, 
en el modelado con el oscilador armónico no se consideran una gran cantidad de 
parámetros que pueden afectar a la descripción del mismo, como lo es la fricción y 
otros. Si consideramos la fricción dentro de nuestro sistema masa resorte, obtenemos 
un movimiento armónico amortiguado cuya solución es más compleja y, de esta forma, 
conforme vamos considerando elementos de nuestros sistemas que pueden o no tener 
impacto en el mismo, nuestro modelo se va complicando y es más difícil de obtener una 
ecuación que lo describa. 


En los últimos años, la implementación de derivadas de orden no entero, es decir 
derivadas de la forma D?, f(x) = En tal que a € Q, ha tenido un gran impacto 
en la descripción de sistemas relacionados con interacciones de largo alcance y sistemas 
con memoria descrita por una potencia [4]. En particular, el cálculo fraccionario ha 
descrito con gran exactitud fenómenos relacionados al amortiguamiento [5]. Como 
consecuencia, el cálculo fraccionario nos ha permitido investigar la elasticidad en 
materiales desde una nueva perspectiva. Para esto se utilizan modelos fraccionarios 
no locales de materiales elásticos, en los cuales el kernel (1 —+t)*% de las integrales 
fraccionarias describen un nuevo comportamiento mecánico del material [6]. En el 
caso de sistemas visco-elásticos, como las arterias, la potencia del kernel de la derivada 
fraccionaria permite pasar de sistemas completamente elásticos a sistemas viscosos [7]. 
Para ilustrar esta última idea, notemos que las arterias suelen ser modeladas mediante 
el uso de Dashpots, los cuales son obtenidos a partir de la ley de Hook y son descritos 


por la siguiente ecuación: 


de(t) 
o(t) = ng 
donde d es el impulso aplicado, 7 es la viscosidad y e es la presión. Los Dashpots describen 
un comportamiento completamente viscoso, por lo cual no se obtiene una descripción 
adecuada de las arterias ya que estos vasos sanguíneos se comportan como resortes. Como 
propuesta ante esta dificultad se introducen los Spring-pot, los cuales se definen como 


o(t) =yD“e(t) 
donde Df, es la derivada fraccionaria de orden «z, tal que n € Q. Este elemento 


interpola un resorte (a: = 0) y un Dashpot (a = 1) [7], permitiéndonos de esta forma, 
describir el comportamiento de las arterias de una forma elegante y sencilla. 
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La versión fraccionaria del operador de Laplace también ha tenido una gran aplicación en 
distintas áreas como mecánica y electrostática, inclusive en el marco de la probabilidad, 
donde actúa como generador de procesos de Levi estable [8]. En el área económica, las 
derivadas de orden no entero se han usado en procesos económicos dinámicos, en el 
cual el orden fraccionario de la derivada modela la memoria del sistema, parámetro que 
suele ser olvidado en la mayoría de los modelos [9]. 


Las derivadas fraccionarias cuentan con propiedades inusuales, como la violación de la 
regla de Leibniz (ley del producto) y de la regla de la cadena (derivada de la composición 
de dos funciones ) [9]. Si bien es cierto que estas propiedades nos ayudan a describir 
características complejas de sistemas dinámicos, dichas propiedades llegan a presentar 
más dificultad en la resolución de ecuaciones diferenciales fraccionarias, por lo cual se 
han utilizado diferentes enfoques para intentar hallar una solución analítica de ecuaciones 
diferenciales fraccionarias como por ejemplo: la aplicación directa de la transformada de 
Laplace, el uso de operadores y el método de análisis con homotopía [10]. 


En este trabajo, se presentan dos métodos de resolución de ecuaciones diferenciales 
parciales fraccionarias. En el primer método, se estudia la transformada de Laplace 
de la derivada fraccionaria de Riemann-Liouville. Posteriormente, proponiendo una 
solución de tipo producto de funciones, se separa la ecuación diferencial parcial 
fraccionaria en dos ecuaciones diferenciales ordinarias fraccionarias [11, 12] las cuales 
son resueltas mediante el uso de la transformada de Laplace [13], de tal forma que la 
solución de la ecuación diferencial parcial fraccionaria original pueda ser expresada 
como producto de las soluciones de las dos ecuaciones diferenciales ordinarias 
fraccionarias. En la segunda metodología, se asume que nuestra ecuación admite 
derivadas sumables, lo cual nos permite integrar en ambos lados de la ecuación 
diferencial parcial fraccionaria y, utilizando las relaciones entre la derivada e integral 
fraccionaria de Riemann-Liouville, se llega a una ecuación Laplace transformable[14]. 
Tomando la transformada de Laplace del resultado previamente obtenido y 
reorganizando los términos se llega a una ecuación de Volterra de segunda especie, 
que garantiza una única solución[15]. 


HISTORIA DEL CÁLCULO FRACCIONARIO EN IR 


El origen del cálculo fraccionario data de 1695, año en el que el marqués de 
UHópital mediante una carta, preguntó a Leibniz cual sería el resultado de E con 
fla) == yn-= 5 A lo que Leibniz respondió que “esto conduciría aparentemente 
a Una paradoja de la cual algún día serán extraídas consecuencias muy útiles”. Con 
estas palabras nació el cálculo fraccionario [16]. A partir de ese momento, muchos 
matemáticos contribuyeron al desarollo del cálculo de orden no entero. En 1730, Euler 
desarrolló trabajos sobre la interpolación entre derivadas de orden entero. En 1812, 
Laplace definió la derivada fraccionaria en base a integrales. Posteriormente en 1819, 
el matemático francés Lacroix discutió por primera vez la derivada fraccionaria en un 


libro de cálculo [17]. A partir de 


Na) =x 
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Lacroix determinó que la n-esima derivada de f(x) = 2” estaría dada por 
ago m! 
T— = 2", m>an. (1) 
dir (m=—n) 
Con la cual, para n = 3 y m = 1 se obtiene el siguiente resultado 
dz 2yz 
dal Jr" 


ha . A _ 2yz . 
De forma similar, si tomamos, n = 3 Y Ha) = e se obtiene 


de? 2/z , 

di? yr 
de lo cual se concluye que ¿e dr de 
dol dol? dx 


En 1822, Fourier fue el siguiente en dar una definición de derivada de orden arbitrario. 
Sin embargo, no se le dio utilidad a esta definición de derivada fraccionaria. El primer 
uso de derivadas de orden fraccionario fue dado en 1823 por Abel [17, 18], en el 
cual aplicó cálculo fraccionario para solucionar la ecuación integral que surge de la 
formulación del problema de la tautócrona. Posteriormente, en 1832 Liouville publicó 
tres memorias sobre el tema [19] junto con algunos trabajos de investigación. En dichos 
trabajos, Liouville definió a las derivadas fraccionarias con el uso de series infinitas; no 
obstante, las derivadas definidas mediante series infinitas estaban limitadas al intervalo 
de convergencia de la serie. Finalmente, mediante el uso de integrales, Liouville llegó a 
construir una mejor definición de derivada fraccionaria, que incluso le permitió tomar 
potencias negativas [17, 19]. 


En 1869 Sonin trabajó inicialmente en la definición de Riemann-Lioville, en un escrito 
llamado “En la diferenciación con índice arbitrario”, empezando con la fórmula integral 
de Cauchy. Letnikov escribió cuatro escritos referentes al tema los cuales tituló “Una 
explicación de los principales conceptos de la teoría de diferenciación de Índices 
arbitrarios” donde dio una extensión de los escritos de Sonin. Sin embargo, no fue si 
no hasta 1884, que el matemático francés Laurent publicó sus escritos de la teoría de 
generalización de operadores partiendo del análisis complejo, dando así las bases para el 
calculo fraccionario [20], y referencias en el mismo. 


PRELIMINARES 


En esta sección se introducirán los objetos necesarios y se desarrollará una notación para 
el desarrollo del trabajo. 


Existen funciones que juegan un papel importante en el desarrollo y entendimiento 


del cálculo fraccionario, a saber, la función Gamma y la de Mitagge-Leffler, las cuales 
juegan un papel importante al momento de buscar soluciones explícitas de ecuaciones 
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diferenciales parciales fraccionarias. Por ende, en las siguientes subsecciones se definirán 
y mencionarán algunas propiedades relevantes de estas funciones. 


Función Gamma 


La función Gamma aparece en problemas relacionados con la normalización de 
funciones de onda de Coulomb o en el cálculo de probabilidades en mecánica 
estadística [2] y juega un rol importante en el cálculo fraccionario. La función Gamma 
está definida como 


TM =/ etéldt, 2€C, Re(z)>0, 
0 


donde Re(z) denota la parte real de un numero complejo z= a + bi, es decir 
Re(z) = a. La restriccion Re(z) > 0 es una condición necesaria para garantizar 
la convergencia de la integral. La función Gamma puede ser interpretada como la 
generalización de un factorial, ya que ésta permite calcular factoriales de números no 
enteros debido a que cumple la siguiente relación: 


T()=(G-1), 


por lo cual P(z) también satisface la siguiente ecuación funcional [21]: 
TG+1) = 21 (2). 


Dada la relación que guarda la función Gamma con los factoriales, se puede generalizar 
la ecuación (1) obtenida por Lacroix de la siguiente manera 


dix” ==, T(m — 1) gon m>n 
dan — Tím=n-1)! y e 


Para mas detalles sobre la funcion Gamma se recomienda ver [2]. 
Función de Mitagge-Leffler 


Una de las funciones mas importante en este trabajo es la función de Mittagge-Leffler, 
ya que a pesar de no ser mencionada en las tablas de la transformada de Laplace, la 
función de Mitagge-Leffler surge en el cálculo de la transformada inversa de Laplace de 
funciones de tipo s“(a+bs%), donde s es el parámetro de la transformada de Laplace 
y a, b son constantes reales [22]. Por otro lado, la función de Mittage-Leffler también 
aparece en la solución de algunos problemas de valores de frontera que involucran 
ecuaciones integrodiferenciales de tipo Volterra con parámetro fraccionario[22, 23]. Está 
definida sobre todo el plano complejo como 


Eaplz) = $ er; a, B,2 E C,Re(a) > 0,Re(B)>0. Y 


n=0 


En base a la ecuación (2) se puede expresar una gran variedad de funciones, por ejemplo: 
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1 
Eo (2) = 21 |] <1, Eri(2)=e*, E2r(z) = cosh(y2), 


iia mas YO iia 


vz 
La función de Mittage-Leffler suele ser vista como una generalización de la función 
exponencial, por lo cual muchas de la soluciones de ecuaciones diferenciales 
fraccionarias son expresadas en términos de la función de Mittage-Leffler, al igual que 
las soluciones de ecuaciones diferenciales ordinarias son expresadas en función de 
senos y cosenos [17]. 


Cálculo fraccionario en R 


A lo largo de la historia han existido varias definiciones para la derivada fraccionaria 
como lo son la derivada de Riemann-Liouville, Hadamard, Erdélyi-Kober, Caputo, Caputo- 
Fabrizio, Marchaud, por mencionar algunas [23, 24]. Sin embargo, para este trabajo 
únicamente nos vamos a enfocar en la definición de Riemann-Liouville, motivados por 
las aplicaciones en diversas áreas como nano electrónica, física y biología [25-33]. 


Definición 1 (Integral Fraccionaria Izquierda de  Riemann-Liouville). Sea 
b(x) € Li(a,b), entonces la integral fraccionaria izquierda de Riemman-Liouville 
(12,9) (x) de orden a. > 0 está definida como 


Maó) = 55 [ At a 6) 


La cual suele ser denotada como (D9)(x). 
Definición 2 (Derivada Fraccionaria Izquierda de  Riemann-Liouville). Sea 


f: [a,b] —+R, entonces la derivada fraccionaria izquierda de Riemman-Liouville 
(De, f)(x) de orden a > 0 está definida como 


(Ds. Ja) = (E) 1 j at 2>0, 


donde n = [a] + 1 y [a] es la parte entera de o. 


Por lo cual, tomando en cuenta la integral fraccionaria de Riemann-Liouville se tiene 

que las derivadas fraccionarias derecha e izquierda de Riemann-Liouville pueden ser 
A _ dyjla a — ¿yla ñ 

reescritas como Di = q7ly" y Dir = lo respectivamente. Notemos que, 

cuando () < a < 1, entonces la derivada fraccionaria izquierda de Riemann-Liouville 


esta dada por la siguiente identidad: 


apa A 1 ELO, o (4) 
(Do, f)(a) = E Es 
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De tal forma que cuando a —> 1, la derivada fraccionaria D2, f(w) se aproxima a la 
derivada clásica L f(x) en x = a.En este trabajo analizaremos ecuaciones diferenciales 
con derivadas fraccionarias cuyo límite cuando y => 1,son derivadas clásicas de orden 
dos. Por lo cual, sin pérdida de generalidad utilizaremos la ecuación (4) para todos los 
cálculos posteriores. Para ilustrar cómo funciona la derivada fraccionaria izquierda de 
Riemann-Liouville vamos a considerar el siguiente ejemplo: 


Ejemplo 1. Evaluemos D*, (1 — a)?, donde B > —1y a >0. 
Considerando la derivada fraccionaria izquierda de Riemann-Liouville enunciada en la 
Definición 2, en particular considerando la ecuación (4) se tiene: 


A d 1 * (tay? 
ado deT(1—a) Sl (2 — ya 


Mediante el cambio de variable z = t — a, tenemos que dz = dt, por lo que 


Realizando otro cambio de variable y = E, tenemos que (x — a)dy = dt, por lo tanto 


nl / (a — ajo te (1 — y) y 


dalla 
e - _—— B(B+1,1-a) 
HO +DO=0 0 q 
= ña E yo 4% 
Finalmente, concluimos que 
cto) E a ay 6) 


Notemos que, sien la ecuación (5)tomamos 8 = O,entonces (x — a)? = (1— a)? =1, 
yas la derivada fraccionaria de una constante k de orden « es 


DU k = == 
Ql 


Como se puede observar, tanto la definicion de derivada como la definición integral 
fraccionaria incluyen integrales. Estas integrales pueden llegar a diverger, por lo que es 
importante considerar el espacio de funciones sobre el cual deseamos trabajar. En este 
artículo vamos a denotar como 12, (L1) a la clase de funciones f representadas por la 
ecuación (3) de función integrable, es decir f = 12, / tal que y € L1(a, b). 
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Teorema 1. Sea a >0. Entonces f E Li+(L1), si y solo si I75*f € AC"(a, b)), 
n= la] +1 y (1759 (a) =0, k=0,n-— 1. Ver la prueba en [23]. 


Observación 1. AC”(a, b] denota la clase de funciones f, tal que f es m— 1 veces 
diferenciable en [a, b] y f(1=D es absolutamente continua en [a, b). 


Cabe mencionar que las funciones pertenecientes a 15, (L1) enfatizan que la 
representación de una función f mediante una integral fraccionaria de orden «a. y la 
existencia de la derivada fraccionaria de orden q: de f son dos cosas distintas [23]. Por 
otro lado, la teoría de Lebesgue nos permite integrar por partes sólo si las funciones 
son absolutamente continuas. Debido a esto, es natural encontrar que la suposición de 
integrabilidad y existencia en casi todo lado de las derivadas fraccionarias (D%, f)(x) 
no es suficiente para nuestro estudio, ya que la integrabilidad y existencia en casi todo 
lado de las derivadas fraccionarias no garantiza la representación de f mediante 
integrales fraccionarias. Por estas razones, para estudiar la reciprocidad entre la derivada 
y la integral fraccionaria es necesario imponer una condición más fuerte sobre f, con lo 
que se da paso a la siguiente definición: 


Definición 3. Sea a: > 0. Una función f € Li(a, b) tiene una derivada integrable Dis f 
siIM7*f(a) € AC”([a, b]), n= [a] + 1. 


De tal forma que la relación entre la derivada y la integral fraccionaria se enuncia en el 
siguiente teorema: 


Teorema 2. Sea a > 0. Entonces la ecuación 


(Di+la+0) (1) =0(1) 
se mantiene para cualquier función integrable p(x). Sin embargo, 


(13,D%. Na) 6) 


sólo se satisface si 
f(x) € L+ (La). 


Si asumimos que f € Li(a,b) tiene una derivada integrable (D2, f)(x) respecto a la 
definición 3, entonces la identidad (6) no se cumple; de hecho, es remplazada por el siguiente 
resultado: 


(12,D9, (a) = (a) + > A Ed a (O 


con n = [a] + 1. Ver la prueba en [23]. 
Proposición 1. (Propiedad de semi-grupo de la derivada fraccionaria). Sea a > 0 


y B>0 tal quen=-1<a<nym-1<B<m y a+B< 1 y sea 
f € Li(a,b) y I"=% € AC"([a, b]). Entonces, se cumple la siguiente relación: 
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(Ds-D?, $) (2) = (DI? 1)(z) — E 


j=1 
Ver prueba en [23, 34]. 


En particular la propiedad del semigrupo es Dear, = De.DE. siempre y cuando 


Pla") =0, j=0,1,---,n—1, 


f € AC”""([a,b]), f” € La(a,b) y n= [8] +1 


Existe la definición de la derivada aproximándose a la derecha, ver [23]. No 
obstante,nosotros sólo nos enfocaremos en la izquierda. 


Transformada de Laplace 


La idea principal de la transformada de Laplace es convertir una ecuación que involucra 
derivadas e integrales en una expresión algebraica [35]. Esto permite resolver ecuaciones 
diferenciales de una forma más sencilla, en especial si se cuenta con condiciones 
iniciales [17]. 


Para ilustrar la relación que guarda la transformada de Laplace con funciones importantes 
del cálculo fraccionario, como son la función Gamma y la función de Mittage-Leffler, 
consideremos las siguientes proposiciones: 


Proposición 2.Sea t > 0, ty = 0ya > —1. Entonces la transformada de Laplace de 
F(t) = t* existe, y está dada por la siguiente ecuación: 


Lt) = T(a +1) 


gal” 


Demostración. Sea t > 0, ty = 0ya > —lentonces la transformada de Laplace de 
F(t) = t* está definida de la siguiente forma 


LAto) =e / te de, (7) 


qe ñ A d A 
utilizando el cambio de variable t = 2, tenemos dt = E por lo cual (7) se reescribe 
de la siguiente forma 


ete=/ Ge. 
0 


Sacando + de la integral y dado que a > —1, se obtiene una función Gamma, con 
z = (+1 de tal forma que 
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pa — a 


Proposición 3. Sea t > 0, tp =0 ys*>a a€R, entonces la transformada de 
Laplace de f(t) =P E,, (at) existe, y esta dada por la siguiente ecuación 


sob 


(8) 


LLE, g(at”)) = 


sea 


Demostración. Sea t > 0, ty =0 ys*>a a € R, entonces la transformada de 
Laplace de f(t) = +4 E,, g(at”) está definida de la siguiente forma 


LG E, ¿(at%)) = / UA E, ¿(at%) de 
0 


Pay 
=/ ena mara t 


n=0 


ar 00 00 
— ¿Qn+B=1 ¿st gg 
Tína + 6) yl 2 


Dada la convergencia de la integral y de la serie podemos intercambiar estas [21]: 


NN Q = q 2 an: ES —$S 
O re Y 


n=0 
[oe] e 
E f an+B=1 
2 Tína + B) E ) 
Z ar T(an+58) (9) 
E 2, T(na +8)  sen+B 


Transformada de Laplace de la Derivada Fraccionaria 
A partir de esta subsección vamos a denotar a las funciones evaluadas en 0, es decir 


F(0), como f(t)li=o. Esto permite involucrar condiciones iniciales al problema. Sin 
embargo, las condiciones iniciales no son discutidas en este artículo. 
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Proposición 4. Sea f una función Laplace transformable. Entonces la transformada de 
Laplace de la derivada fraccionaria de orden a: de Riemann-Lioville centrada en a+ = 0 
está dada por 


L(DG-$()) == LAO) - DL (0) limo: eds 
Demostración. consideremos la derivada fraccionaria de Riemann-Liouville 


a _ 4, 1 AO 
(Des J(x) da ( +] E-ga% 


l-a 


la cual puede ser reescrita como 


DF = 90 


(=a +1) 


donde G(t) = fit 7) ef(mar. Tomando la transformada de laplace de (11) 
obtenemos 


AO = AC) 
= 5 e(G(0)- Gto 


elf torio 


= Ao) 


= Tan A 7 


Noa (=a +1) 
2 £  Fro+1 _ E ltlimo_ 
—T=a+1) sia Or T(—a +1) 


cor (y [tm leo 


Cambiando e: por —a + 1 y usando la ecuación (3) obtenemos que 


LIDGF (0) = LF) - DO limo . n 


Proposición 5. Sea n(x) una función Laplace transformable. Entonces la transformada 
de Laplace de la derivada fraccionaria de orden 2a: de Riemann-Liouville centrada en 
ar =0 es 


LLDBEn(e)) =s“L(m(t) — s*D ¿LU n(t)limo — Dj n(t)limo. 


Demostración. Sea n(x) una función Laplace transformable. Entonces por la 
Proposición 4 se obtiene 
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L(Dim(t) =s*LEDgn(t)) — DL Dgntt)limo: 


Simplificando, 


LíDBEn(e)) = se (s*Lin(t)) — DL la) limo) — Dl 
= LA n(t)) =D ¿Etico — DL ft)lio: 


+ 
== 
= 

11 

o 


Procediendo de igual forma que en la demostración de la Proposición 4 y utilizando las 
propiedades de convolución, se pueden mostrar las siguientes relaciones: 


* £ (72, u(o,1)) =s"U(z,s), 

ko L (E no O 
* £L (Ue tiuta.0) = E OS dy, 
e ala) ) 00) 


Donde F;(s) = L4f,(t)). 
Espacio de Zemanian 


Se conoce al espacio de Zemanian como el espacio en el cual viven las transformadas de 
Laplace. Por lo cual, para obtener la transformada inversa de Laplace se necesita tomar 
en cuenta su forma generalizada en dicho espacio. Por lo que tenemos las siguientes 
relaciones presentadas en el trabajo de Ferreira y Vieira [14] 


Ly 030. 1 SO ods = (DAL q.) (8), 


—ir 


a+ir 
lm+30 / ¿”+6 gslt—to) q — gl (+ to). 
O 


ir 
donde ¿=0,1,...,n n € No, d esla función delta de Dirac, y su convergencia es en D' 
Ecuaciones integrales 


Las ecuaciones integrales son aquellas que relacionan una función incógnita con una 
integral, de tal manera que la función incógnita se encuentre bajo el signo integral. 
Análogamente a las ecuaciones diferenciales, las ecuaciones integrales pueden ser de 
orden lineal y no lineal. En el caso de ecuaciones integrales lineales podemos encontrar las 
de tipo Fredholm y las de tipo Volterra [15]. 


1. Ecuaciones de tipo Fredholm: Son aquellas ecuaciones integrales con límites de 
integración fijos. Éstas pueden ser clasificadas en dos grupos: 
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Primera especie: La función incógnita sólo está bajo el signo integral: 
1 
Ha)= [ Kanéwdy 052 ysl 
0 


donde p(x) es la función incógnita, K(x, y), f(w) son funciones conocidas y K(x, y) es 
continua en el cuadrado O < x, y < 1 [36]. 


Segunda especie: La función incógnita se encuentra bajo el signo integral y fuera 
de él. 


A / El éllda, 05% YyEL hell 


Si f(x) = 0, se conoce como ecuación homogénea, mientras quesi f(w) 4 O seconoce 
como ecuación completa o no homogénea. 


2. Ecuaciones de tipo Volterra: Son un caso particular de las ecuaciones de tipo Fredholm 
en el cual se satisface que K(x, y) > 0 si y > zx. Por lo cual, de la misma forma que con 
las ecuaciones de tipo Fredholm, podemos clasificar las ecuaciones de tipo Volterra en dos 
grupos 

Primera especie: La función incógnita sólo está bajo el signo integral. 


Ha) = / Klan Es YEn 


Segunda especie: La función incógnita aparece tanto bajo el signo integral como 
fuera de él. 


T 
Ha) =0(2)-A [Ka oda 05% ysl AEC. 
0 
Observación 2. A la función K(x, y) sele conoce como el Kernel o núcleo de la ecuación. 
A pesar de la complejidad de las ecuaciones integrales, el método de descomposición de 
Adomian [15, 37] es efectivo al momento de buscar una solución de ecuaciones integrales. 


Dicho método consiste en descomponer la función incógnita en una serie infinita de la 
siguiente forma: 


donde Po(1) = f(x) y f(x) son determinados a partir de la siguiente fórmula recursiva 


Alisa / EA 


Utilizando el método de Adomian se puede encontrar una solución para ecuaciones 
integrales de Volterra de segundo tipo. En particular, dado un kernel K(x,t) = (1 —t)07! 
se puede enunciar el siguiente teorema. 
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Teorema 3. [15] Sea f € L1(la, b]), a > 0 y A € C. Entonces la ecuación integral 


A 


ula) = 8(0) + a / (0H u(pdt, zea. 


Tiene una única solución dada por 


i= e / 0 EVaMa — 9 F (dt. 


CÁLCULO FRACCIONARIO PARA ECUACIONES DIFERENCIALES 


En este trabajo se estudian dos métodos distintos para encontrar soluciones explícitas 
de ecuaciones diferenciales parciales fraccionarias, las cuales fueron propuestas por 
grupos de trabajo interdisciplinario. Para esto, se considera el siguiente conjunto de tres 
ecuaciones diferenciales parciales, que son parte central de los resultados obtenidos por 
dichos grupos 


Dil, t) +Cu(x,t) = Deule, £), ecuación tipo evolución, - (13) 
(0) 
Dizu(z, t) +Culz,t) = D”.u(z, t), ecuación tipo onda, (14) 
0 
Dizu(z, D+ Di¿U(z, t) = Di. u(z, t), ecuación tipo mixta, (15) 
10) 


de tal manera que cuando a > 1y B > 1, las ecuaciones (13), (14), (15) se comportan 
de forma clásica como se muestra a continuación: 


E t) +Culz,t) = Cule,t) CER, 


dx 

de d 

qa, t) +Culzx,t) = qe, t), CER, 
d? d d 
qa (e, t) + Co Az, t) = qa, t), CER. 


Solución por producto de funciones 


La metodología usada en esta sección consiste en reducir una ecuación diferencial 
parcial fraccionaria a dos ecuaciones diferenciales ordinarias fraccionarias que puedan 
ser resueltas mediante la transformada de Laplace, de tal forma que la solución del 
problema inicial sea el producto de éstas. Para esto, vamos a considerar la derivada de 
Riemann-Liouville de orden 24 centrada en a+ = 0, esto es, (D22). La construcción 
de las soluciones por esta metodología puede ser resumidas en los siguientes pasos: 
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Paso 1: Proponer una solución de tipo producto de funciones u(x,t)= u(x)u(t) tal 
que u(x), u(t) € C[O, b], y separar la parte espacial de la temporal. 


Paso 2: Resolver las ecuaciones diferenciales fraccionarias homogéneas mediante la 
transformada de Laplace. 


Paso 3 : Expresar la solución del problema original como un producto de funciones. 


Para ilustrar dicho método, consideremos la siguiente ecuación diferencial parcial 
fraccionaria: 


Défu(x,t) + Cu(z,t) = Dé, u(z,t), (16) 


Comencemosconelpaso1.Consideremos u(x,t) = u(x)u(t) talque u(x), u(t) € C(0, b] 
(a, B) € [0,1] x [0,1], (x,t) € Q = [0,b] x [0,b] y b < oo. Reemplazando u(x, t) 
en la ecuación (16) y realizando algunos cálculos directos se tiene que 
1 


ER 2222 u(a)] = — |D% u 
0 Pétulz) +Culz)] = oz [Dit] 


Para que los dos lados de la ecuación (17) sean iguales, estos deben ser igual a una 


constante, por lo cual se tiene que 


Du(a) + (C — Au(x) =0, (18) 
Dé, u(t) — Mít) =0. (19) 


(17) 


Continuando con el paso 2, buscamos soluciones para las ecuaciones (18) y (19) con 
A € R. Estas soluciones pueden ser obtenidas mediante las siguientes proposiciones 


Proposición 6. Sea u: [0,b] +R, AE R,y0 < a < 1, entonces la ecuación diferencial 
fraccionaria ordinaria 


Di+ula) + Pula) =A, PER, 20) 


tiene como solución 
ula) = A Evan (PIS) + DA llaot Esa EP.  Q1) 


Demostración. Sea u : [0,0] +R, AER, y 0<a< 1, por lo cual tomando la 
transformada de Laplace de (20) se tiene 


LiDg+u(e)) + PLfu(z)) ==, 


= SLlula)) Dt fallo + PL) = a 


Ordenando la ecuación (22) y mediante algunos cálculos directos se tiene la siguiente 

solución 

2 A y Der olaaa 03) 
s(so + P) (so + P) 


Lfu(z)j 
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Finalmente, con ayuda de la ecuación (8) con B=a>+1, B=a0 para el primer y 
segundo término, se toma la transformada inversa de Laplace de la ecuación (23). Así 
obtenemos 


ue) = A Enjara (Pa%) + Dl) leon Evo (Pat). m 


Proposición 7. Sea u : [0,0] +R, AER, y 0<a< l, entonces la ecuación 
diferencial fraccionaria ordinaria 


AD%u(w) + Bu) =A, A,BERN (O), (24) 


tiene como solución 


A B (1-20 a Ha 
u(u)= Ga Bao (-q2) + ae MD aleoz? as (Ge ) 


== == B 104 
+ De Vu(a)lamo2” “Ba (2 ) . (25) 


Demostración. Sea u: [0, b] >Ry0Q <a < 1. Entonces la ecuación (24) puede 
ser reescrita como 


Ditu(o) +w*%u(x) =c, (26) 


donde w? = ES y c= A . Tomando la transformada de laplace de (26): 


L[D¡tu(o)) + w*L1fu(u)) = Lfc), 
= s*Lfu(a)) - “DL ulo)lao — Dilo) lao +? Lfu(x)) = - 


0+ 


La cual puede ser escrita de forma reducida como 
Lula) A - A+ Lula) => 27) 
(1-20 


donde A, = DA lalo y As = Do+ ula)lemo. Ordenando la ecuación 
(27) y mediante algunos cálculos directos se obtiene que 


Cc + Ás ; Aso (28) 
($20 + 2) * ($204?) * (s20 4?) 


Llu(e)) => 


Finalmente, con ayuda de la ecuación (8) con B=2a>+1,P=2a0, 6 =0 para el 
primer, segundo y tercer término respectivamente, se toma la transformada de Laplace 
del lado derecho de la ecuación (28). Así obtenemos 


ula) = a Exogosr (2190) + Aga? Eso ga (0) y Aya Ego (tao). 


O equivalentemente: 
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A B (120 ” Bo, 
u(z) = Ga Erazo (42) + PRA Ñ u(a)la=ot? Esa za (qe ) 
ik _ B 
+ De u(2)loo2*" Exa (qe) s a 


Para concluir con el paso 3, usando las Proposiciones 6 con A= 0y P = Ay 7 con 
A=1,B=(C-—A) y A=0 obtenemos que la solución de la ecuación diferencial 
parcial fraccionaria (16) está dada por 


u(z,t) = (Dla)lao Bao zo (=(C - AJa?o) 
HD (alamo Eno ( (0 — Ma?) 


x (Delia E (-A19)) 


De tal forma que se puede enunciar el siguiente teorema: 


Teorema4.Seau(x,t) = u(x)ju(t) talqueu(x)u(t) e Cl0,b), (a, B) € [0,1] x [0,1] 


,(x,t) € Q =[0,b] x [0,b], b < oo, entonces la ecuación diferencial parcial fraccionaria 
Déu(x,t) + Cu(z,t) = Dé, u(e, t), 
tiene como solución 
u(z,t) = [AO TN ((C- AJa%) 
4D Aa lao2 zoo (=(C— yaa) 


x (Dto Ea. (-A19)) l 


donde C € R, (a, B) = (0, 1] Xx [0, 1), (x,t) € 0 = [0, b] AS (0, bl, db < 00 


De manera análoga, usando apropiadamente las Proposiciones 6 y 7, se demuestra el 
siguiente teorema: 


Teorema 5. Sea u(x,t) = u(av)u(t) tal que u(x) u(t) € C[0,b] (a, 8) e [0,1] x (0, 1), 
(x,t) e€ Q=[0,b] x [0,b], b < oo, entonces la ecuación diferencial parcial fraccionaria 


Di+u(zx,t) + Culz,t) = Df, u(z, da 
tiene como solución 
u(z,t) = (DL) la=02 Esa (=(C — a")) 


x (DL altlimot Ea (=A19)), 


donde C' e R, (a, 6) € [0,1] x [0,1], (x,t) e Q= (0, b] x [0, b], b< oo. 
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Solución por integración directa 


La metodología usada en esta sección consiste en integrar directamente la ecuación 
diferencial parcial fraccionaria, de tal forma que, tomando en cuenta la relación entre 
la integral y la derivada fraccionaria, se encuentra una ecuación cuyos términos 
son Laplace transformables. Para esto, vamos a considerar (a, 8) € [0, 1] x[0, 1] 
(x,t) € Q = [2o, Xo] x [to, To], Zo,to > 0, Xo, To <oo, tal que u(x,t) admite 
derivadas fraccionarias integrables. La construcción de las soluciones por esta 
metodología puede ser resumida en los siguientes pasos: 


Paso 1: Se integra directamente y se toma en cuenta la relación entre la integral y la 
derivada fraccionaria (Teorema 2). 


Paso 2: Se aplica la transformada de Laplace para llegar a una ecuación integral de tipo 
Volterra de segunda especie. 


Paso 3: Se resuelve la ecuación integral y se simplifica. 
Paso 4: Se toma la transformada inversa de Laplace. 


Para ilustrar dicho método consideremos la siguiente ecuación diferencial parcial 
fraccionaria: 


De u(a,t) + Culx,t) = DP ,u(z,t). (29) 


Comencemos con el paso 1. Consideremos la ecuación (29), donde CER, 
(a, B) € [0,1] x [0,1], (2,t) € Q = [2o, Xo] x [to, To] , To,to > 0, Xo, To < oo, y 

u(x,t) admite derivadas fraccionarias sumables ED y Dias Aplicando el operador 
integral > a los dos lados de la ecuación (29) y, tomando en cuenta la relación entre 
la integral y la derivada fraccionaria enunciada en el Teorema 2, se tiene que 


1 
TA u(zo, 1) +CIhPu(z,t) = 0 e t). 


alo : B e ; 
De forma similar, aplicando el operador 1;,+ a los dos lados a la ecuación anterior y 
utilizando el teorema de Fubini, se toma en cuenta la relación entre la integral y derivada 
fraccionaria enunciada en el Teorema 2 de tal forma que se obtiene 


zp)" ,, TY a a 
Pala) EC De ula) - A DSulao,t) + CIS 12 (e, t) 
tt) 
= LA ufa,e E mo Li u(z, to) 


la cual es equivalente a 
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= Ea u(x,t) — 0 Lo go(2), (30) 


donde fo(t) = Il, %u(zo,t), fu() = De,.ulzo,t) y gol) = TC u(s, to), que son 
las condiciones de Cauchy. 


Continuando con el paso 2, se toma la transformada de Laplace respecto a £ de la 
ecuación (30) y utilizando sus propiedades de convolución se obtiene que 


— ío)* =p je -B ó 
sU(x, s) == = sE (s) (2 o sPEo(s) + a E (x _ y “Uy, 5)dy 
Ñ mes 1) / E (31) 


Multiplicando la ecuación (31) por sé y reorganizando sus términos obtenemos la 
siguiente ecuación de tipo Volterra: 


En == Zg) ( = zo)" 6 —tos pa+l 
U(z, s) TA+ yal )+ T(a) Fo(s) — e La Y (2) 62) 
sB— E 
+ E ll (1 — y) U (y, s)dy, 


Procediendo con el paso 3, se toma en cuenta que la ecuación (32) en virtud del Teorema 
3, tiene como única solución 


U(.0) = E o) + EE o) te rz(a) 


4 (578 - C) 1KE 0 Y" Eotta+1 ((s* => C)u => y)”) 


TO 


(Po) + E to) — rta) de 


(33) 


Dada la convergencia de las integrales y las series en la ecuación (33) podemos 
intercambiar éstas y reescribirla de la siguiente forma: 


20 
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U(z, t) = (x se a ((s* > Cia > 2o)2**) Fo(s) 
+ (2 — 20) Ea+1,0+1 ((s* - C)x — 20)0**) F1(s) 


00 
Y e A pla) 
n=0 


Para concluir con el paso 4, se requiere cancelar la transformada de Laplace, por 
consiguiente tomamos en cuenta su forma distribucional y con ayuda del teorema 
multinomial y mediante algunos cálculos directos, podemos dejar la solución expresada 
en términos de la función de Mittage-Leffler 


ula,t) = (020) Easra (Dijo — O) — 20)%**) fo(s) 
+ (020) “Barros (WD, — Ola — 20)**1) f(s) 


= [towns (Dé, ON 19") 5 talud 


ZO 


Donde go(x) = TI” Pu(z, to), fo(t) = To u(zo,t), f,(t) = D%,u(xo,t)y $ es la 
función delta de Dirac. 


De tal forma que se puede enunciar el siguiente teorema 


Teorema 6. Sea CER, (a, 6) € [0,1] x [0,1], (x,t) € Q = [2o, Xo] x [to, To] 
Zo,to >0, Xo, To < oo, tal que u(x,t) admite derivadas fraccionarias integrables 
on 7 Fs + Entonces la ecuación diferencial parcial fraccionaria 


Du(x,t) + Culz,t) = Df u(z,t), 


xo+ 


tiene como solución 


ula,t) = (020) Enrra (Df. — ON 20)%+) fo(s) 
+ (0 20) Eagrarr (Dé, — C)lo — 20)%*) (5) 


= [Eneas ((D%, — Ole") 50 1o)anl0dda 


0 


donde go(1) = Luz, to) folt) = I¿“u(zo,t), f(t) = Disulzo,t) y 8 es la 
función delta de Dirác. A 
Análogamente se puede demostrar el siguiente teorema 


Teorema7.SeaC € R(a, B) € [0,1] x [0,1] (x,t) € Q = [zo, Xo] x [to, To), To, to > 
0, Xo, To < oo, tal que u(x, t) admite derivadas fraccionarias integrables D£_, y De 
Entonces la ecuación diferencial parcial fraccionaria 


Di +u(z,t) + Culx,t) = DF +u(x,t), 
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tiene como solución 
u(e,t)= (2) Ena ((;,. - Ola — 20)”) folt) 
= [00 Ena (0%, ON 9)") 5 toa, 


TO 


donde go(w) =1'; a Bula, to) fot) = en “u(zo, t) y $ es la función delta de Dirac. 


CONCLUSIONES 


Si bien el uso de la transformada de Laplace es un gran método para resolver 
ecuaciones diferenciales, éste puede involucrar cálculos extremadamente difíciles. A 
pesar de que una gran cantidad de transformadas inversas de Laplace de funciones tipo 

(a +bsf%) pueden ser expresadas mediante la función de Mittage-Leffler, el cálculo 
de transformadas inversas de Laplace está limitado por el orden de «+, ya que al ser éste 
variable, es difícil encontrar una factorización adecuada mediante fracciones parciales. 
Dicho esto, se puede observar que proponiendo una solución de tipo producto 
de funciones podemos llegar a encontrar de una forma más sencilla la solución de 
una ecuación diferencial parcial fraccionaria de segundo orden, reduciéndola a dos 
ecuaciones diferenciales ordinarias fraccionarias homogéneas. A pesar de ello, la 
propiedad de semi-grupo de la derivada fraccionaria de Riemann-Liouville nos dice que 
para derivadas de orden 2a nuestra solución solo es válida para a € (3, 1). 


Por otro lado, se puede observar que la relación entre la integral y derivada fraccionaria 
puede llevarnos a encontrar ecuaciones de Volterra de segunda especie. Si bien es 
cierto que la importancia del orden de integración y derivación depende del espacio 
en donde se trabaje, el uso de funciones Ly nos permite llegar a ecuaciones de Volterra 
de segunda especie, asegurando así una única solución. Hay que destacar que la 
existencia de la ecuación de tipo Volterra depende de que todas las derivadas respecto 
a la misma variable tengan el mismo orden, es decir no se puede tener derivadas de 
orden Den y Dis dentro de una misma ecuación diferencial parcial fraccionaria si se 
desea buscar una ¿olución por el método de integración directa. La ventaja de usar el 
método de integración directa es que permite extenderlo a múltiples dimensiones de 
una manera mucho más simple que el método de la solución por producto. 


La particularidad de estos métodos es que están limitados a coeficientes constantes 
por lo cual el siguiente paso a seguir sería trabajar con ecuaciones diferenciales 
fraccionarias, parciales y ordinarias con coeficientes variables. Para esto, se propone en 
un futuro definir de una manera conveniente la ley de la cadena y expandir al método 
de Frobenius al caso fraccionario, ya que muchas de las ecuaciones diferenciales que 
describen fenómenos físicos son resueltas por estos métodos. En particular, el método 
de Frobenius nos ayuda a resolver ecuaciones diferenciales cuya solución son funciones 
especiales, como los polinomios de Hermite o Laguerre [2], de tal forma que extendiendo 
el método de Frobenius al caso fraccionario podamos calcular funciones especiales de 
orden no entero. 
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Finalmente, este trabajo se limitó a buscar soluciones generales a las ecuaciones 
diferenciales parciales fraccionarias propuestas, por lo que no se profundizó sobre 
problemas de valores iniciales y problemas de valores de frontera específicos sobre los 
resultados obtenidos. Sin embargo, las soluciones obtenidas consideran las condiciones 
iniciales y condiciones de frontera como constantes, de tal forma que los grupos de 
trabajo interdisciplinario puedan usarlas con mayor libertad. 
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